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We study the one-dimensional SU(4) exchange model under magnetic fields, which is the simplest 
effective Hamiltonian in order to investigate the quantum fiuctuations concerned with the orbital 
degrees of freedom in coupled spin-orbit systems. The Bethe ansatz approaches and numerical 
calculations using the density matrix renormalization group method are employed. The main concern 
of the paper is how the system changes from the SU(4) to the SU(2) symmetric limit as the magnetic 
field is increased. For this model the conformal field theory predicts an usual behavior: there is 
a jump of the critical exponents just before the SU(2) limit. For a finite-size system, however, 
the orbital-orbital correlation functions approach continuously to the SU(2) limit after interesting 
crossover phenomena. The crossover takes place in the magnetization range of 1/3 ~ 1/2 for the 
system with 72 sites studied in this paper. 
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I. INTRODUCTION 

In recent years increasing attention has been paid to 
the physics related to the orbital degrees of freedom and 
its quantum fluctuations in the strongly correlated elec- 
tron systems in addition to the charge and/or spin de- 
grees of freedom. In the conventional treatments for the 
coupled spin and orbital systems, mean-field-type ap- 
proaches have been used to determine an ordered orbital 
structure and the effects of the quantum fluctuations are 
usually neglected |l]j|]. Recent experiments suggest, how- 
ever, in some transition metal and rare earth compounds, 
the orbital degeneracies and quantum fluctuations of spin 
and orbital degrees of freedom play an important role in 
their unusual properties, which may require more sophis- 
ticated treatments than the mean- field approximations. 

For examples, CeBg is known as a typical dense Kondo 
system {Tk = 1 ^ 2 K) with the Fg quartet ground 
state, and is well known for its unique magnetic field- 
temperature phase diagram Since the quadrupolar 
ordering temperature (Tq ^ 3.4 K) is of the same or- 
der as the spin Neel ordering temperature {T^ 2.3 K), 
the interplay between spin and orbital quantum fluctu- 
ations may be important and thus we need to consider 
the effects of quantum fluctuations more seriously beyond 
the mean-field theories. A pure LiV204, which has a 
fee normal-spinel structure and does not show any static 
magnetic order, is found to be the first compound which 
shows the heavy fermion like behavior among the tran- 
sition metal oxides Q. But its origin of the large effec- 
tive mass may be different from the usual heavy-fermion 
systems which have nearly localized /- and itinerant con- 
duction electrons. The highly frustrated lattice structure 
and nearly degenerate t2g orbital at each V'^ '^+ ion are 
believed to be the origin of the unusual properties, such 
as disordered ground state and heavy-fermion like behav- 
iors. 

Generally speaking, the quantum fluctuations are less 
important in higher dimensions, and the classical pic- 
tures, i.e. spin and orbital ordered structures in the case 



of the couples spin and orbital models, may be applied 
for many cases in three dimensions. However, these ex- 
periments suggest the possibilities that the interplay be- 
tween spin and orbital fluctuations destroys the classical 
ordered structures and the ground state may become a 
disordered liquid state. 

In order to study the physics in which quantum fluc- 
tuations concerned with the orbital degrees of freedom 
are important, the SU(4) exchange Hamiltonian (|l|) has 
been considered as a prototypical model and studied ex- 
tensively by several groups This SU(4) exchange 
model is derived from the highest symmetric limit of the 
two band Hubbard model with the strong correlation 
{U/t 3> I) and the vanishing Hund coupling (J = 0) 
at quarter filling. The model is given by. 
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where {i,j) stands for a nearest neighbor pair and K is 
given by 2t^/{/ > 0. It consists of usual S- and pseudo 
T- spins, which describe the spin and orbital degrees of 
freedom, respectively. It is important to note that the 
both spins have the Heisenberg type isotropic interac- 
tions and the total symmetry of this model is the SU(4) 
group, which is higher than the apparent SU(2)(g)SU(2) 
group. 

In one dimension, overall properties of the ground 
state and excitations of this most symmetric SU(4) ex- 
change model are well understood by the previous studies 
[H,0-|To|. In terms of the spin-ladder models, Eq.(|l|) corre- 
sponds to the two-leg AF Heisenberg ladder with 4-body 
interaction producing frustrations. A similar model with- 
out frustration (the 4-body interaction with the opposite 
sign) is studied by using the green's function Monte Carlo 
method . It is noteworthy that the present model be- 
longs to the Bethe solvable class with the general SU(A^) 
symmetry and also belongs to the quantum critical sys- 
tems with strong quantum fiuctuations. 
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In actual systems, the highest SU(4) symmetry may 
be destroyed into various lower symmetries due to some 
anisotropies and different physics emerges depending on 
the way how the symmetry is lowered. In many phys- 
ically relevant systems, the spin rotational invariance is 
approximately hold, while the SU(2) symmetry of the or- 
bital part is easily destroyed. For example, it becomes 
U(l) by a single- ion anisotropy and a finite Hund cou- 
pling makes the exchange interaction of the orbital part 
Ising type with Z2 symmetry Another type of sym- 
metry lowering is the models with SU(2)(8)SU(2) sym- 
metry which is given by Eq.(|l]) with arbitrary constants 
replaced for 1/2. This model in one dimension is in- 
vestigated both numerically and analytically by several 
authors 

In this paper we study the one-dimensional SU(4) ex- 
change model under magnetic fields as another example 
with a lower symmetry. It is naturally expected that the 
quantum fiuctuations in one dimension are more promi- 
nent than those in two or three dimension. Therefore one 
dimensional model may shed right on the role of quantum 
fiuctuations concerned with the spin and orbital degrees 
of freedom. It should be also mentioned that this model 
is equivalent to the case with a splitting between the 
two orbitals. Under magnetic fields, the highest SU(4) 
symmetry is broken by the Zeeman splitting and the to- 
tal symmetry of the model becomes SU(2)(g)U(l), where 
the SU(2) symmetry represents the rotational invariance 
of the orbital degrees of freedom and the U(l) symme- 
try describes the spin degrees of freedom under uniaxial 
magnetic fields. In sufficiently strong magnetic fields, it 
is evident that the all spins are fully polarized and our 
model reduces to the SU(2) symmetric antiferromagnetic 
(AF) Heisenberg model with only pseudo T- spin degrees 
of freedom. 

Therefore, it is interesting to investigate how the sys- 
tem changes from the SU(4) to the SU(2) symmetric 
limit as we increase the magnetic field. We study this 
transition by calculating the ground-state energy and el- 
ementary excitations and also by looking at the behav- 
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iors of the orbital-orbital correlation functions. For this 
purpose, we have used numerical calculations with the 
density matrix renormalizing group (DMRG) method as 
well as analytic approaches based on the Bethe ansatz 
and the conformal field theory (CFT). Since these meth- 
ods do not introduce any uncontrolled approximations, 
they are suitable for investigating the effects of quantum 
fiuctuations beyond the level of various mean-field theo- 
ries. 



II. GROUND STATE AND EXCITATIONS 
A. Bethe ansatz equations and symmetry properties 

We consider the one dimensional SU(4) exchange 
model with the periodic boundary conditions (PBC). 
First we calculate the thermodynamic ground-state en- 
ergy under magnetic fields by using the Bethe ansatz 
methods for finite and infinite size systems. In the ac- 
tual calculations, we solve the Bethe ansatz equations 
with specified values of the magnetization m^, which is 
defined by 2Sfgf/N. In this notation, rriz — corre- 
sponds to the SU(4) symmetric limit with no magnetic 
field and = 1 corresponds to the SU(2) symmetric 
limit where all S- spins order in parallel and only orbital 
degrees of freedom are left. The ground-state energies 
for the SU(4) and SU(2) symmetric limit are obtained 
analytically to be 1 - 7r/4 - 3 log 2/2 ~ -0.825 118 9 and 
1 - 21og2 ~ -0.386 294 4, respectively §. 

Let us suppose that iVi, iV2, N3, and iV4 are the 
numbers of four species {N — Ni + N2 + + N4 
and N is the total number of sites) with the quantum 
numbers (5^T^) = (T,T), (T,i), (i,T), and (1,1), re- 
spectively. These numbers specify the value of S^^^ — 
{N,+N2-N3-N4)/2 andT^^t = {Ni-N2 + N3-Ni)/2. 
We also define Mi ^ N2 + N3 + N^, M2 = N3 + N^, and 
M3 = N4. The Bethe ansatz equations for the SU(4) 
exchange model with the quantum numbers A^i, N2, N3, 
and A'4 are given as follows [pj: 
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where K{a) = 2tan~-'^(2a) ± tt, 9{x) = — 2tan~-'^a;, and 
the suffix and k are numbered from 1 to Mi, M2 and 
M3, respectively. For a given set of Ja,Jf}, and we 
can calculate the parameters ai , (3j and 7^ , so-called ra- 
pidities, by solving the Mi + M2 + coupled non- linear 
equations (|^). The quantum numbers J a, Jp and J-y 
are integers (half integers) depending on the value of the 
M1+M2, Ml + M2 + M3, and M2 + M3 to be odd (even), 
respectively, and their distributions specify a state. The 
ground-state energy per site (sat), momentum (P), and 
rriz for the specified number of species, TVi, A^2,-^3j and 
N4, are given as follows: 

2 Ml Ml ^ 

1=1 i=i » 

Ml / Ml M2 Ms \ 

^ = E^M = ^ E^^+E-^/^.+E-^- h (4) 

ruz = ni + n2 - ns - n^, (5) 

where we take K as energy units hereafter, putting K — 1 
in Eq.(|l|), and define Tij = Ni/N. In order to avoid 
complications coming from multi-fold degeneracies 
we only consider the system with N — -in sites (n is an 
integer) in the following. 

Under magnetic fields, the four- fold degenerate bands 
in terms of the original Fermion model, are split into two 
doubly degenerate bands by the Zeeman effect and these 
two bands are labeled by spin up and down, respectively 
(see Figj^). Since S- spins are gradually polarized with 
the increase of magnetic fields and pseudo T- spins are 
not affected directly, it is reasonable to suppose that the 
ground state in magnetic fields is represented by a irre- 
ducible representation with a finite S^^^ and zero T^^j, 
which is given by the Young's diagram as illustrated in 
Fig.^. For TTiz = and 1, the irreducible representations 
of the ground state are given by the SU(4) and SU(2) 
singlet, respectively. For < rriz < 1, the upper and 
lower half rectangles represent the local SU(2) symmetry 
of the orbital degrees of freedom for the up and down 
spin bands split under magnetic fields, respectively. 




FIG. 1. Schematic band structure of the SU(4) exchange 
model under magnetic fields, ef denotes the Fermi energy. 
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FIG. 2. The Young's diagrams for the ground states un- 
der magnetic fields. Below each diagram, magnetization and 
its total symmetry are given. 



B. Ground-state energies under magnetic fields 

We have numerically solved the coupled equations (||) 
under the conditions, A^i — N2 and N-^ = iV4, and by 
using Eq. (|^) and (|^) the ground-state energies with spec- 
ified magnetization rri^ are calculated for the system size 
N — 100, 200, and 400, respectively. In order to esti- 
mate the bulk limit, we have extrapolated to the infinite 
system size limit by fitting these data with the function 
ejv = eoo+ci/N + C2/N'^. Bulk ground-state energies per 
site as a function of are obtained as shown in Fig.^by 
the solid line. The ground-state energy for the SU(4) and 
SU(2) symmetric limit are calculated to be —0.825 119 
and —0.386 2944, which are perfectly consistent with the 
analytic results. 

-0.35 I • • • • 1 




FIG. 3. Thermodynamic ground-state energy as a func- 
tion of rriz . The solid line represents the ground state-energy 
obtained by the extrapolation to infinite system size. Broken 
lines stand for the asymptotic behaviors around the SU(4) 
and SU(2) symmetric limits. 

In order to consider the thermodynamic properties of 
the SU(4) exchange model and to analytically obtain the 
asymptotic forms around the SU(4) and SU(2) symmetric 
limits, we take the infinite system-size limit in the Bethe 
ansatz equations under the conditions of A^i = N2 
and = By introducing distribution functions per 
unit length for the continuous rapidities a, /3, and 7, we 
obtain the simultaneous integral equations as follows: 
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where G„(a;) = n/7r(l+n^x^) and i?i, i?2, and i?3 are the 
distribution functions per unit length for the rapidities a, 
/?, and 7, respectively. The infinite cut-ofF parameters ex- 
cept for those of /3 in the integral equations preserve the 
relations A^i = N2 and A^3 = A^4 and the cut-ofF B de- 
termines the difference between A^i and N3 and specifies 
the magnetization rriz through the relation below. The 
ground-state energy per site (soo), the momentum P and 
the magnetization are given by 

£00 = 1 — 27r / G2{a)Ri{a)da, (7) 
P = N [ K{a)Ri{a)da, (8) 

J-B 

m,=2 [ R2{(3)d(i. (9) 

J\13\>B 

By solving the coupled integral equations (^) numeri- 
cally and using the equations and (||), we may again 
obtain the thermodynamic ground-state energy as a func- 
tion of rriz- Since there exist non-systematic errors in 
the procedures of discretizations for the integral equa- 
tions, it is hard to calculate the ground-state energies for 
the SU(4) and SU(2) symmetric limits in better accuracy 
than 0.1% relative errors to the exact values. We have 
found that, for the calculation of the bulk ground state 
energies, Bethe ansatz method for finite-size systems with 
the aid of the extrapolation is more suitable than that for 
the infinite-size system. Since the SU(4) and SU(2) sym- 
metric points correspond to the limits with the cut-off 
parameter B — 00 and 0, respectively, we get analyti- 
cally the asymptotic expressions of the thermodynamic 
ground-state energies (Eoo) around these limits within the 
lowest-order calculations as follows: 

^00-1-5-^ + -^^^^, (for^.^O), (10) 
4 2 (log tan 212^) 

- 1 - 21og2 + log2 (m^ - 1) , (form^~l). (11) 
We also show these functions in Fig.0 by the broken lines. 



C. Excitations and incommensurate soft modes 

In order to obtain excited states by the Bethe ansatz 
method for a finite-size system, we follow the standard 
approach to change the distributions of integer or half- 
integer quantum numbers Jq, Jg, and Jj and calculate 
the elementary excitations with a specified TOz. A re- 
moval, from the symmetric and continuous distributions 
of the quantum numbers for the ground state, creates 
a hole and an addition creates a particle. By creating 
particle-hole pair excitations corresponding to the three 
kinds of quantum numbers Ja, JfS, and Jj , we can discuss 
three elementary massless excitations in our model. The 
calculated excitation spectra as a function of momentum 
q are shown by the solid lines in Fig|^(0)~(3) for several 
values of to^. We have defined the thermodynamic limit 
by the extrapolation of finite-size data in the same way 
as for the ground-state energies and defined the momen- 
tum of the elementary excitations by the relative one to 
the ground state. 

Figure ^ shows how the excitation spectra change be- 
tween the SU(4) and SU(2) symmetric limits. In Fig.§(0), 
we see that three excitations have the same velocity tt/2, 
which is consistent with the exact result obtained by 
Sutherland With the increase of magnetizations, 

the amplitudes and softening wave vectors of the two 
small branches of the excitations are decreased continu- 
ously and disappear at the SU(2) limit, while the largest 
branch approaches to the SU(2) elementary excitation. 
For niz = 1, the model reduces to the SU(2) AF Heisen- 
berg model of the pseudo T- spins and its elementary 
excitations have been obtained to be tt sin q, which are 
nothing but the des Cloizeaux and Pearson modes. It is 
important to note that for accurate calculations of the 
smallest elementary excitations for larger m^, we need 
calculations for larger system sizes because those excita- 
tions correspond to the rearrangements of the quantum 
number J-y and its number M3 = N(l — mz)/4: is de- 
creased with the increase of ruz ■ In other words by using 
the Fermion model, this difficulty comes from the fact 
that the relative weight of the down spin band to the up 
spin one is decreased with the increase of to^ for a fixed 
system size. 

The wave numbers which correspond to the incommen- 
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surate soft modes in the elementary excitations (Fig.^) 
are specified by the nesting vectors of the Fermi points 
in Fig.|l| as follows: 2kF,i — (1 ~ "^z) /2, 4fci?_| = 
7r(l — ?7iz), and 2kp^i = tt (1 + TOj) /2. Since the rela- 
tion 4:kp^i = 4fci?_| (mod. 2it) is always satisfied, the 
4fci?_| mode is not independent. We have found that 
the excitation spectrum with the softening at q — 2fci?j 
(2fci?.|) correspond to the kink-type excitations of the or- 
bital pseudospin T in the down (up) spin band and the 
spectrum with the softening at q — Akp^i corresponds to 
the kink-type S- spin excitations, so-called spinous. 

It may be worth to mention that the shift of the 
momentum of the singularity in the structure factor 
is actually observed for the SU(2) spin case. Ben- 
der et.al. have studied the copper benzoate, 
Cu(C6D5COO)2-3D20, which is known to be the typ- 
ical spin 1/2 one-dimensional antiferromagnet, by the 
neutron scattering experiment and have found the exis- 
tence of the field-induced soft modes at the wave vector 
q — n±2TTM{H)/gfiB, which is consistent with the Bethe 
ansatz solutions. 




t t 31/2 f It/2 

2kF,i 4kF,\ 2kf i 



ZKfj 

FIG. 4. Dispersion relations in the infinite-size system 
limit for several values of m^. The vertical and horizontal 
axis represent the energy and momentum, respectively. Be- 
low each figure, we label the incommensurate soft modes in 
terms of the original fermion model. 



III. CROSSOVER PHENOMENA OBSERVED IN 
ORBITAL-ORBITAL CORRELATION 
FUNCTIONS 

A. Correlation functions in real space 

In this section, we investigate the orbital-orbital corre- 
lation functions under magnetic fields by the numerical 
calculations using the DMRG method Since the 



DMRG method is more suitable for open boundary con- 
ditions (OBC) than PEC, we employ the OBC in order 
to carry out the calculations for large system sizes with 
sufficient accuracy. We consider the 72-site system and 
obtain the correlation functions up to the distance of 36 
sites. As usual we calculate the wave function of the 
system with a specified S^^f instead of dealing with the 
model under magnetic fields. To eliminate the oscilla- 
tions coming from the open boundaries, we define the 
averaged orbital-orbital correlation functions as follows: 

(^^^A-,)^ jEm+feT^4fe+,> • (12) 

fe=0 

We have confirmed that after this averaging procedure 
the i dependences of the correlation functions are almost 
removed. The obtained correlation functions are shown 
in Fig.|for = 0,1/3,1/2,5/6, and 1. 

For general models with the SU(Af) symmetry, the 
asymptotic behavior of the correlation functions have 
been obtained in a simple form by Affleck jlOj. When 
we apply his results to our model with the SU(4) or 
SU(2) symmetric limit, which corresponds to — or 
1, respectively, the orbital-orbital correlation functions 
behave asymptotically as follows: 

/T-T- \ A ^ cos(27ri/iV-K^o) 1 \ 

^^ u+j) - ^0 1^ -p ' 

where Aq and Lp^ are non-universal constants depend- 
ing on each model, while the critical exponents a — 
2(1 — 1/iV) and (3 = 2 are universal constants. For the 
SU(4) symmetric limit, we try to fit the data from j = 18 
to 32 with this asymptotic form by the least mean square 
method. We do not use the data for j > 32, because these 
sites are so close to the boundary that they are affected 
by the oscillations from the open boundary even after the 
averaging procedures. From the fitting we get the critical 
exponents to be a = 1.52 ± 0.08, and /3 = 2.02 ± 0.08. 
These exponents are consistent with the analytic results 
by Affleck. For the other limit with the SU(2) symmetry, 
we similarly estimate the critical exponents for the oscil- 
lating term with wave vector q = n and non-oscillating 
term to be 1.04±0.01 and 2.5±0.6, respectively, neglect- 
ing the logarithmic corrections. These results are again 
consistent with the exact values 1 and 2. In this fitting 
analysis, we have used the data from j = 18 to 25 in order 
to eliminate the boundary effect, which is considered to 
be more noticeable than that for the SU(4) limit because 
the critical exponent of the SU(2) limit is smaller than 
that of SU(4). 

For the intermediate magnetizations between the 
SU(4) and SU(2) symmetric limits, we can not get mean- 
ingful results concerned with the exponents from a sim- 
ilar fitting as the SU(4) or SU(2) limit, since the cor- 
relation functions have more complicated structure as is 
shown in Fig.^. However, the correlation functions in 
real space show interesting behaviors as the SU(2) sym- 
metric limit is approached. First of all, let us take a look 
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at Fig.|| paying attention to their periodicities and how 
they approach to the SU(2) hmit. For a small value of 
niz, we see clear periodic oscillations with the wave vec- 
tors 2kp^i and 2kp_^ (see FigJ^(a) and (b)), which cor- 
respond to the incommensurate soft modes of the exci- 
tation spectrum discussed in the previous section. As 
we increase magnetic fields, the correlation functions ap- 
proach to the SU(2) one (Fig.||(e)) in such a way as the 
number of periodic structures observed within the 36- 
site length is decreased and at the same time the SU(2)- 
like structure with the wave vector tt becomes noticeable 
(Fig.^(c) and (d)). These behaviors are consistent with 
the fact in the elementary excitations, that 2kF,i ^^'^ 
ikp^i soft modes vanish and the 2fcF,t soft modes contin- 
uously evolve to the des Cloizeaux and Pearson modes 
at TT with the increase of magnetic fields. We have also 
found that for a large value of rriz, the correlation func- 
tions are well-illustrated by the following pictures where 
the oscillations with the wave vector tt are modulated by 
the slowly-varying 2kp^i^. 

With the increase of m^, the period of the 2kp^i oscil- 
lations becomes longer and the region, where the SU(2)- 
like behaviors with wave vector tt are dominant, becomes 
more extended, and in the end for = 1, the correlation 
functions become the true SU(2) ones. In this way, the 
bigger the value of is, the more distant becomes the 
asymptotic region, which is governed by the envelope of 
2kp^l oscillations, and the more prominent are the SU(2)- 
like behaviors within the envelope. Generally speaking, 
as long as we consider a finite-size system, the calculated 
correlation functions do not reach the asymptotic region 
for a larger value of than some critical magnetization, 
which is an increasing function of the system size, and 
the crossover phenomena from the SU(4) to the SU(2) 
symmetric limit take place. Beyond this characteristic 
magnetization the SU(2)-like behaviors with the slowly- 
varying envelope described by the 2kp^i oscillations are 
a better picture for the calculated correlation functions 
in a finite length than the asymptotic behaviors of the 
2kp,i and 2kp -^ oscillations. 
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FIG. 5. Orbital-orbital correlation functions (TlTf^j) as a 
function of j for the 72-site system with = 0, 1/3, 1/2, 5/6 
and 1. In figures (a) and (e), we plot the fitting results by the 
symbols +. 



B. Correlation functions in momentum space and 
asymptotic behaviors 



Since we want to minutely observe the transition from 
the SU(4) to the SU(2) symmetric limit with the increase 
of the magnetic field, next we focus on the asymptotic 
behavior of the correlation functions and concretely esti- 
mate the crossover region between these two limits for a 
finite-size system with 72-site length. For this purpose, 
we discuss the anomalies of the structure factors by cal- 
culating the fourier transformation which is defined by 



N/2 
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j = -N/2+l 



(13) 
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In this definition, we have periodically arranged the nu- 
merically calculated correlation functions up to N /2-sAe 
length to form the ring with the N sites, which are less 
affected by boundary effects, and defined the structure 
factors for the ring. The obtained structure factors for 
the 72-site system by changing from to 1 by 1/6 are 
shown in Figj^. In the figure we see clear two singular 
cusp structures at g = 2kp^i and 2/cF,t- These domi- 
nant cusps correspond to the field-induced soft modes in 
the dispersion relations obtained in the previous section, 
and may be interpreted as the orbital excitations in the 
down and up spin bands as illustrated in Fig.^. It is 
important to note that the q = 4A:f.| (= 4fci?j (mod. 
2tt)) mode does not appear in the orbital-orbital cor- 
relation functions. This is because the excitation with 
q = '^kp^i is considered to be an excitation with the spin 
currents. In fact the '^kp^i (Akp^-^) momentum transfer 
comes from such excitations as we move two particles 
from the left Fermi point of the down (up) spin band to 
the right one. This situation is completely analogous to 
the case of the one-dimensional Hubbard model, where 
the both 2k p and Akp oscillating parts appear in the 
density-density correlation function, while the Akp part 
is absent in the spin-spin correlation functions because 
the dominant oscillating term with q = Akp originates 
from the holon excitations [1161. 
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FIG. 6. Fourier transform of the orbital-orbital correla- 
tion functions for = 0, 1/6, 1/3, 1/2, 2/3, 5/6, and 1. We 
label the wave numbers 2fci?.| and 2fci?.| by the down arrows 
(I) for each m^. T^ijv) for = 1 (which is out of the plot 
range) is equal to 1.486 . 

According to the conformal field theory the crit- 
ical exponents ai, a2, and which correspond to the 
oscillating terms with q = 2kp,i, Akp^i, and 2kp,^, re- 
spectively, are obtained as follows: 



ai = as 



2 



a, a2 = = 4(a - 1), 



where ^ is defined by ^(/3 ~ B) and ^(/3) is obtained by 
the following integral equation: 



m = 1+1 dp' £,{(3') 

' -B 



\x\ 



1 + el^ 



c(/3-'3') 



We calculate these critical exponents as a function of , 
where niz is given by Eq.(|^), which only depends on B. 
The obtained results are shown in Table |. In the table, 
we represent the limit to unity from < 1 by 1-0. The 
critical exponent a is equal to 1.5 at rriz = and 1-0 and 
continuously changes between < < 1 with larger 
values about 1.5 ^ 1.6 and jumps to 1 at the SU(2) sym- 
metric point with niz = 1. In Fig.^, the left and right 
cusp structures correspond to the modes with the critical 
exponent ai and 0^3, respectively. Though the conformal 
field theory predicts the same critical exponent for the 
both cusps, the calculated structure factors seem to be 
quite different, especially for large m^. Apparently the 
left side cusps become less dominant, while the right side 
cusps are more prominent as the increase of . 

In order to consider possible reasons, we again em- 
ploy the schematic band structure of the original fermion 
model shown in Figj^. As the increase of applied mag- 
netic field, it is evident that the relative weight o f the 
down spin band is diminished as mentioned in Sec.IIC. 



It is the reason that the left hand side cusp structures 
are suppressed and their amplitudes become smaller in 
contrast to the predictions by the CFT. By a closer look 



at the right side cusp structures, it is found that until 
the value of less than about 1/3, the cusp shows a 
little weaker anomaly than that for niz — 0, which is 
consistent with the prediction by the CFT, while for the 
region of > 1/2 the cusps grow rapidly and seem to 
continuously develop into the SU(2) AF cusp structure 
with the critical exponent 1. 

These characteristic features of the fourier transformed 
orbital-orbital correlation functions for the system with 
72 sites are explained by the crossover phenomena due 
to the finiteness of the system and do not contradict the 
CFT predictions. It is reasonable to consider that in the 
region > 1/2 the relative weight of the down spin 
band is so small and the 2kp^i oscillations have such a 
long period that the calculated correlation functions do 
not reach the asymptotic region within the 72-site length 
and show the properties of the single up spin band which 
contribute relatively large weight, that is the SU(2)-like 
behavior in the envelope of the 2kp^^ oscillation. There- 
fore the cusp structures seem to continuously evolves into 
the behavior of the SU(2) AF Heisenberg model. On the 
other hand for rriz < 1/3 both upper and lower split 
bands have significant weight and 2fci?,x and 2kp^-\ are 
not so small that the calculated correlation functions are 
in the asymptotic region and they show the asymptotic 
behaviors of the SU(4) exchange model under the mag- 
netic field predicted by the CFT. This picture is also 
consistent with the behaviors of the real-space correla- 
tion functions in which the SU(2)-like behaviors become 
noticeable around = 1/2. 

To summarize, for small the orbital-orbital correla- 
tion functions show the oscillations with the wave vectors 
2kp,^ and 2fci?,|, which give a consistent picture with the 
prediction by the CFT. For the understanding of the re- 
gion with large niz, however, we find it useful and appro- 
priate to rewrite these two oscillating term as cos{2kp^p) 
and cos(2fci?jj) = cos(7rj) x cos{2kp^p). The second 
equation, which comes from the relation 2kp^i+2kp^i; = tt 
originating from the quarter filling, gives for the correla- 
tion functions of large niz a simple picture in which the 
SU(2)-like oscillations with the wave vector tt are modu- 
lated by the slowly-varying 2kp^i oscillations. Following 
this picture, we may understand straightforwardly how 
a finite-size system continuously makes a transition from 
the SU(4) to the SU(2) symmetric limit as we increase 
magnetic fields. In particular, the continuous growth of 
the cusp structure into the SU(2) limit shown in Fig.|^ 
does not contradict the prediction by the CFT: a jump 
of the critical exponent between niz = 1-0 and 1. 

TABLE I. The critical exponent a for the magnetizations 
niz from to 1 by 1/6. 



TTlz 





1/6 


1/3 


1/2 


2/3 


5/6 


1-0 


1 


a 


1.500 


1.562 


1.581 


1.585 


1.573 


1.543 


1.500 


1.000 
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IV. CONCLUSIONS AND DISCUSSIONS 

In summary, we have investigated the SU(4) exchange 
model, which describes the strong quantum fluctuation 
limit of the two band Hubbard model, under magnetic 
fields and observed how the system transforms from the 
SU(4) symmetric limit with zero magnetization to the 
fully polarized SU(2) limit by using the Bethe ansatz 
method and DMRG calculations of the orbital-orbital 
correlation functions. First of all, we have kept the 
trace of the ground-state energy as a function of 
and also obtained the analytic forms around the SU(4) 
and SU(2) symmetric limits. In order to study low en- 
ergy properties of this model, we calculated three inde- 
pendent elementary excitations by the standard Bethe 
ansatz technique and found the field-induced incommen- 
surate soft modes. Through the DMRG calculations of 
the orbital-orbital correlation functions under magnetic 
fields, we have found that the 2kp^^ and 2kp -^ modes, 
which carry the orbital pseudospin currents, give rise to 
dominant cusp structures, while the 4:kF,i mode, which 
carries the spin current, does not appear in the structure 
factors. We also conclude that for a finite- length sys- 
tem the crossover phenomena take place in the different 
way predicted by the CFT, and after that the correla- 
tion functions are properly described by the SU(2)-like 
behavior with the slowly- varying envelope of the 2kp,i 
oscillations. As ruz approaches to unity, the period of 
envelope becomes longer and longer and thus the cal- 
culated finite-length correlation functions seem to con- 
tinuously change into the SU(2) ones for a finite size. 
Corresponding to this fact, the cusp structures at the 
the 2k F,^ observed in the structure factors develops into 
the SU(2) ones continuously without showing any jumps 
of the critical exponents. By comparing the calculated 
results and the CFT results, we estimate the crossover 
region between the SU(4) and the SU(2) symmetric lim- 
its to be 1/3 '--^ 1/2 when we consider the relatively short 
length system of 72 sites. It should be noted that from 
the symmetry of S- and T- spin degrees of freedom, we 
expect similar results for the spin-spin correlation func- 
tions under uniaxial pressures or crystal fields produced 
by single-ion anisotropies. Under the anisotropics, one 
type of orbital of the two is preferred and the four-fold 
degenerate bands split into the two doubly degenerate 
bands depending on the = ±1/2. 

Lastly, we would like to briefly comment on the rel- 
evance of the apparent difference in the crossover phe- 
nomena between the calculated correlation functions and 
the correlation exponents predicted by the CFT. Let us 
take the example of rUz = 5/6 shown in Fig||(d). The 
2kF,i = 7r/12 oscillation has the period of 24 sites, and 
thus it is reasonable to assume that it require a few hun- 
dreds sites to reach the asymptotic region. From the CFT 
analysis, we can get the informations about the mathe- 
matically exact asymptotic behaviors for the infinite sys- 



tem size. In the real materials, however, perfect spin sys- 
tem cannot be realized and even for very pure systems 
their intrinsic length, where the exact one dimensional- 
ity are kept without any structural defects or impurities, 
would be some hundred atomic lengths at most. There- 
fore the results discussed in the present paper serve a 
simple example where the crossover phenomena in real 
systems may be different from the one predicted by the 
CFT. In this sense, it is necessary to pay some atten- 
tion for the application of the CFT results to realistic 
finite-length systems. 
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